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ABSTRACT 

The colour-magnitude relation (CMR) of cluster elliptical galaxies has been widely 
used to constrain their star formation histories (SFHs) and to discriminate between 
the monolithic collapse and merger paradigms of elliptical galaxy formation. We use a 
ACDM hierarchical merger model of galaxy formation to investigate the existence and 
redshift evolution of the elliptical galaxy CMR in the merger paradigm. We show that 
the star formation history (SFH) of cluster ellipticals predicted by the model is quasi- 
monolithic, with only ^10 percent of the total stellar mass forming after z ~ 1. The 
quasi-monolithic SFH results in a predicted CMR that agrees well with its observed 
counterpart in the redshift range < z < 1.27. We use our analysis to argue that the 
elliptical- only CMR can be used to constrain the SFHs of present-day cluster ellipticals 
only if we believe a priori in the monolithic collapse model. It is not a meaningful tool 
for constraining the SFH in the merger paradigm, since a progressively larger fraction 
of the progenitor set of present-day cluster ellipticals is contained in late-type star 
forming systems at higher redshift, which cannot be ignored when deriving the SFHs. 
Hence, the elliptical-only CMR is not a useful discriminant between the two competing 
theories of elliptical galaxy evolution. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
formation - galaxies: fundamental parameters 



1 INTRODUCTION 

Evolution with redshift of fundamental physical re- 
lations can provide robust constraints on the epoch 
of formation and subsequent evolution of early-type 
galaxies. The apparently universal relationship between 
colour and luminosity of elliptical galaxies, usually 
referred to as the colour-magnitude relation ( CMR) , 
was first established by ISandagc & Vishvanatharij il977l) . 
although the correlation between these two quanti- 
ties had been demonstrated before (e.g. [B aum| [l959l : 
Ide VaucouleursllQGltlMcClure fc van den Bergbll968h . The 
observed CMR has been widely used as a discrim- 
inant between the two competing theories of early- 
type galaxy evolu t ion, the monolithic colla pse model 
(e.g. iLarsonl Il974l: iKodama fc Arimotol Il997^ an d the 
hierarchical merger model (e.g. iKauffmann et al.l Il993t 
Somorville & Primack 'l999'; 'Co le et al.ll200d: iHatton et alJ 
200.1: .Khochfar fc B urkert 2Q03,) . 

A comprehensive study of the CMR, using photo- 
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metric data based on CCD observations of the nearby 
Virgo and Coma cluster s, was first undertaken by 
iBower. Lucev fc Ellis! (Il99l hereafter BLE92). Their results 
showed a remarkably small scatter about the mean rela- 
tion. Their interpretation of the results, in the context of 
the monolithic collapse model, was to attribute the slope of 
the CMR to a variation in mean metallicity with luminos- 
ity and to attribute the small scatter to a small age dis- 
persion between galaxies of the same size. They concluded 
that the epoch of formation of elliptical galaxies should be 
at z > 2. Subsequent studies of the CMR extended the 
BLE92 results to intermediate redshifts (0 < 2 < 1) and 
showed that there was no detectable evolut i on of the slope 
and scatt er with time (e.g. lE Ujj_gt_glJ|jjgjygt^^ et alJ 
Il998; G l adders et"all ll998l : Xan Dokkum et al.ll200ol) . The 
results from these studies were interpreted as confirmation 
of a high-redshift formation epoch of cluster ellipticals fol- 
lowed by passive evolution to present day. 

Subsequent studies indicated that the key charac- 
teristics of the CMR (slope and scatter) that were de- 
rived by these authors needed some modification. Ellip- 
tical galaxies commonly display radial colour gradients 
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fe.g. Ide Vaucouleurj|l96ll: ISand aee &i VishvanatharJIlQTsI: 
iPranx et alJll989HPeletier et alJll99d) . being optically red- 
der at their cores than at the outskirts. Galaxy colours in 
the majority of these CMR studies were derived using fixed 
apertures which, given that a galaxy's intrinsic size may 
vary, meant sampling different portions of different galax- 
ies. Accounting for the effect of fixed aperture photometry 
is essential in the context of this study, since our model lacks 
spatial information on the scale of galaxies. Producing a cen- 
tral CMR is therefore not possible within the model and we 
must compare our results to aperture-corrected photometry. 

Numerous authors have attempted to quantify the fixed- 
aperture bias that may result due to the presence of colour 
gradients. An efficient way to control this effect is to mea- 
sure the colour insi de an aperture wh i ch scales with the size 
of the galaxy (e.g. iBower et aDll998l: iTerlevich erai]l200lt 
^codcggio 20(^, such tha t one samples an ide ntical fraction 
of the light in each galaxv. feower et alJ il998l) compared the 
CMR slope observed by BLE92 with the slope derived after 
replacing a fixed aperture with the parameter _D„ - the size 
of the galaxy within which the mean surface brightness is 
19.80 mag arcsec"^. They estimated that colour gradients 
accounted for roughly 30 percent of the slope i.e. the mag- 
nitude of the slope in the Dv CMR was roughly two-thirds 
of that derived using fixed-apertures in BLE92. Crucially, 
their study indicated that the CMR maintains a significant 
slope even after correcting for colour gradients. A series of 
other studies have also studied the effect of removing the 
fixed-apertu re bias . Shallow er CMR slopes have been re- 
ported by ; Pru g niel fc Si mien (1996,) who used colours de- 
rived within the effective radius (radius which contains half 
the ga laxy's light. Re) and by iFioc fc Rocca-Vohnerang e 
^9£ ) who used total magnitudes and colours . IScodeggic 



I 20011) suggested that re-computing colours over for the 
BLE92 galaxies causes the apparent slope to decrease from 
-0.082 ± 0.008 to -0.016 ± 0.018, a value which is statis- 
tically consistent with a zero slope. In addition, the scat- 
ter increases from 0.035 to 0.136, due to the large intrinsic 
scatter in the colour gradients. However, the uncertainty in 
colour measurements within R^ can be significantly larger 
than within Dv, beca use the sur f ace b rightness within Dv 
is higher. At 2a the IScodeggiol i200lfl data is consistent 
with slopes between -1-0.02 and -0.52. Modifying the BLE92 
slope for colour gradients, using the 30 percent correction 
factor derived by Bowerjyt al. (1998), gives —0.054, so that 
within the errors there is agreement between the various 
studies. 

Over the last decade, there has been steadily accumu- 
lating evidence for morphological evolution amongst clus- 
ter galaxies, which suggests that formation mechanisms of 
cluster ellipticals are at least not uniquely monolithic. Al- 
though approximately 80 percent of galaxies in the cores of 
prese nt day clusters have early-type morphologies (Dressier. 
Il980l) . a higher fraction of spiral g alaxies have been reported 
in clusters at 0.3 < z < 0.8 (e.g. iButcher fc Oemleiill984l: 
[ Dressier et"ai]|l997l: ICouch et aLlll998l : Ivan Dokkum et alJ 
I2OOOI) . along with increased rates of merger and interaction 
events (e.g. ,Co uch ct a l. 1998; van Dokkum ct al. 1999). 
iKauffmann et ami996l) suggested that only approximately 
one-third of e arly-type galaxies in the Canada-France Red- 
shift Survey dSchade et all Il995^ were fully formed and 
evolving passivelv. iFranceschin^t al.l lll998f) found a re- 



markable absence of early- types galaxies at z > 1.3 in a 
K-band selected sample in the Hubble Deep Field. These re- 
sults strongly suggest that early-type galaxies in nearby and 
distant cluste rs may have been formed from late-type pro - 
genitors (e.g. iButcher fc Oemleilll98 iDressler et al]ll997^ 
and highlight the possible if not essential role of merger and 
interaction events in the formation of early-type galaxies. 
In particular, if the merger paradigm is correct then late- 
type progenitors of the present-day cluster ellipticals must 
be included in any method (e.g. the CMR) employed to de- 
termine their SFHs. Excluding these late-type progenitors 
would produce a distorted view of their formation histories 
(progenit or bias), a point first sugge sted and explored in 
detail bv lvan Dokkum fc Fran:3 (120011) . 

Given the accumulating evidence for formation of 
early-type galaxies from star-forming progenitors at fairly 
recent epochs, a number of authors have success- 
fully reconciled the observed CMR with galaxy merg- 
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Bower et alJll998l: 


Shiova fc Bekki 




van Dokkum et alJ 


Il998l). Apart from 


Kauffmann fc Charlod Jl998l). these studies have not in- 



volved a fully realistic semi-analytical galaxy formation 
model which incorporates the important effects of galaxy 
merging on the chemo-photometric evolution of galaxies. 
One of our aims is to extend these studies by applying a 
ACDM hierarchical merger model to study the CMR from 
low to high redshift. 

We begin our study by discussing the comparative ef- 
fects of age and metallicity in determining the model (U — V) 
CMR at present day and tracing the bulk SFHs of cluster 
ellipticals as a function of redshift. We then explore the pre- 
dicted evolution of the CMR to high redshifts, compare with 
existing observational evidence and, in particular, quantify 
the effect of progenitor bias. Using our analysis of progeni- 
tor bias, we present arguments to show that the commonly 
used elliptical-only CM R, even when i t is derived using equal 
light fractions ( c.f. .Bower et alJll998l : ITerlevich eral]|200ll : 
IScodeggidf200ll) . can only be used to constrain the SFHs 
of cluster ellipticals if we believe a priori in a monolithic 
collapse model. It is not a meaningful method of constrain- 
ing the SFH in the hierarchical merger picture. Hence it is 
also not a useful discriminant between the two competing 
theories of galaxy evolution. 



2 MODEL PARAMETERS THAT AFFECT 
THE PRESENT-DAY CMR 

The model we use in this study is GALICS, which com- 
bines large scale N-body simulations with simple analytical 
recipes for the dynamical evolution of baryons w i thin d ark 
matter haloes. We direct readers to lHatton et al.l ll200 3h for 
specifics regarding the model. There are certain key param- 
eters in the model that affect the age and metallicity of the 
model galaxies and thus have an impact on the slope, scatter 
and absolute colour of the predicted CMR. A discussion of 
these model parameters is necessary, not only to elucidate 
their effect on the CMR, but also because the actual setup 
we use in this study is slightly diff erent from the reference 
model given in lHatton et al.l i2003i) . The setup has been al- 
tered, firstly to make some corrections to the metallicity of 
fresh gas injected into DM haloes, and secondly to bring the 
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predicted metallicities of the model galaxies in agreement 
with current observational evidence. Table^summarises the 
changes in the characteristics of the CMR due to variations 
in these parameters. In the subsequent sections we present 
an explanation of the parameters and the values used in this 
study. 



2.1 Baseline metallicity 

The reference model in iHatton et al.l i2003l) adds pristine 
i.e. metal-free gas to DM haloes when they are identified. 
However, the haloes are not identified until they achieve a 
threshold mass of 10^^ Mq. In reality, early population II 
stars would already have polluted the ISM in the time that 
it takes for such halo identifications to take place. Hence, 
the gas in the haloes should not be pris tine but slightly 
pollu ted. Chemical enrichment models (e.g. IPevriendt et alJ 
Il999r) suggest that this pollution should be of the order of 
O.IZq. Hence, we use this value as a baseline metallicity for 
fresh gas injected into DM haloes in the model. The baseline 
metallicity has a negligible impact on the slope and scatter 
of the CMR but slightly affects the absolute colour of the 
cluster sample as it changes the average metallicities of the 
model galaxies. 

2.2 Threshold black hole mass and IMF 

Another parameter that affects the metal input into the 
ISM, and therefore the average metallicity of the stellar pop- 
ulation, is the threshold mass at which a star becomes a 
black hole (BH). This is still poorly unde rstood, but esti- 
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mates suggest masses around 50 ± WMq iTsuiimoto et all 
119971) . based on a combination of local stellar [0/Fe] abun- 
dances and chemical enrichment analysis. We use a threshold 
black hole mass of 60Mq in this study. 

Since massive stars make a significant contribution to 
the metal enrichment of the inter-stellar medium (ISM), the 
proportion of massive stars and hence the IMF affects the 
mean metallicities of the model galax :ies. In this study we 
use a Kennicut IMF (|K ennicutlj|l98j) , whic h was also used 
in the fiducial model of iHatton et all J2003h . Both the BH 
threshold and IMF increase the dynamical metal enrichment 
of the ISM. This changes not only the mean metallicity of 
the galaxies, and hence their absolute colour, but also makes 
the slope of the CMR steeper. This is because more massive 
galaxies, which have deeper potential wells, retain gas and 
therefore metals more effectively, leading to higher enrich- 
ment of the ISM and stellar populations that are born from 
it. However, less massive galaxies tend to lose their gas con- 
tent anyway, so that a larger injection of metals into their 
ISM does not have a big impact on the metallicity of their 
stellar populations. As a result of this differential behaviour 
the slope of the CMR becomes steeper. 



3 PROPERTIES OF PRESENT-DAY CLUSTER 
ELLIPTICALS 

In this section, we explore the present-day CMR predicted 
by our model. Figure Q presents the predicted CMR in our 
model for cluster ellipticals at z — 0. Also shown is a linear 
least-squares fit to the points (dotted line) and a progressive 




Figure 1. The model colour magnitude relation at 2 = 0. TOP: 
CMR sequence with a linear least-squares fit (dotted line) and a 
progressive one-sigma fit to the sample. BOTTOM: Cluster ellip- 
ticals split into their individual clusters. 



one-sigma fit to the sample, with the error bars indicating 
the local spread of points about the best-fit relation. We 
select present-day cluster ellipticals by identifying elliptical 
galaxies in dark matter (DM) haloes with masses of W^'^Mq 
and above. Also shown in Figure is the CMR sequence 
with galaxies coded according to their mean metallicities 
and ages. The bottom panel splits the model ellipticals into 
their individual clusters. The model slope is derived in all 
cases using a linear least-squares fi t. The scatter is ca lculated 
using Tukey's Biweight statistic ijBeers et al.lll99(il i. which 
has commonly been used by observers in CMR studies. Table 
1 compa res the model CMR with those derived by BLE92 
and Bower et al.l il99ft) within Dv and by^Scodcaizio (20oJ) 
using the effective radius, Re. 

The predicte d model slo p e is c onsistent with both the 
value reported by IScodeggicl i200lfl and the BLE92 value, 
after correcting for colour gradients using the 30 percent cor- 
rection derived in Bower ct al. ( 1998). The predicted scatter 
is smaller than that derived bv lScodeggi o ( 2001) but r oughl y 
1.5 times larger than that reported bv lBower et alj il99di . 
We note, however, that the scatter in the model galaxies it- 
self varies from cluster to cluster, so that the intra-cluster 
scatter may be different from the global value across all clus- 
ters. In Figure|5|we split our sample of model ellipticals into 
their respective clusters and plot the intra-cluster scatter as 
a fraction of the global scatter in the sample. We find, for 
example, that the cluster with the largest number of ellip- 
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Table 1. Variations in CMR slope and scatter with baseline metallicity, threshold BH mass 
and IMF. The slope is derived from least squares fits (over the magnitude range My = — 19 to 
My = —23) and the scatter is calculated using Tukey's bi-weight statistic. The values in bold 
indicate the parameters used for this study. ^BHT = Black Hole Threshold mass. 

Baseline metallicity = 





BHT = 45M0 


BHT = 6OM0 


rirl i = IzUlVlQ 


Kennicutt IMF 


-0.037 ±0.007 
0.072 


-0.049 ± 0.007 
0.079 


-0.050 ±0.009 
0.094 


Scalo IMF 


-0.034 ±0.006 
0.057 


-0.033 ± 0.009 
0.070 


-0.047 ±0.010 
0.10 




Baseline metallicity = O.IZq 
BHTt = 45M0 BHT = 60A/q 


BHT = 120Mq 


Kennicutt IMF 


-0.036 ±0.009 
0.075 


-0.047 ±0.010 
0.082 


-0.052 ± 0.010 
0.12 


Scalo IMF 


-0.032 ±0.007 
0.061 


-0.034 ±0.010 
0.078 


-0.045 ± 0.010 
0.10 



Table 2. Comparison between the characteristics of our model CMR at ^ = w i th BL E92, 
corrected for col our gradients using the 30 percent correction given in lBower et alj il99St) and 
IScodeggiol 120011) who used the effective radius {Re) of galaxies to derive colours. For the model 
CMR the slope is derived from least squares fits (over the magnitude range My = —19 to 
My = —23) and the scatter is calculated using Tukey's bi-weight statistic. 

Source Slope Scatter (mag) Colour 

Bower et al. (1992. 1998'! (Coma) -0.054 ± 0.007 0.049 U-V (within D„) 

Scodeggio (2001') (Coma) -0.016 ± 0.018 0.136 U-V (within i?e) 

This study -0.047 ±0.010 0.082 U-V (total magnitudes) 

ticals has a lower scatter than the global value, although 
there are other large clusters which exhibit a scatter above 
the global value. We note that our model galaxy sample is 
an ensemble of galaxy sets from different clusters, whereas 
the observations usually refer to only one cluster. There may 
be additional issues contributing to the discrepancy between 
the model and observed scatter - for example, there appear 
to be strong radial colour gradients in cluster population s at 
low redshift (e.g. de Propris et al. 2004, MarEonincr c t alJ 
l200ll : lEilingson et al.l200ll) . such that bluer objects reside in 
the outer parts of clusters. The scatter of the observed CMR 
would therefore depend on the maximum cluster-centric ra- 
dius sampled in the observations. In addition, since the ob- 
servations are not derived from total colours, it is possible 
that colour gradients are correlated with deviations from the 
true total-colour CMR in such a way that aperture colours 
appear to have smaller scatter. For these reasons, we do not 
find the discrepancy between the model and observed scat- 
ters particularly compelling. 

We find that the model predicts a significant correla- 
tion between colour and luminosity. An important question 
is how the model CMR is generated at present day. Clearly, 
in a hierarchical merger scenario age is expected to play a 
part in generating any such sequence. It is therefore crucial 
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Figure 2. Intra-cluster scatter as a fraction of the total scatter 
plotted against the elliptical occupancy of each cluster. 



to disentangle the effects of age and metallicity and deter- 
mine how much of the correlation is generated by a variation 
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Figure 3. Variation of mass-weighted average ages and metallic- 
ities with absolute V-band luminosity and baryonic mass. 
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Figure 4. Age-metallicity parameter space in model cluster ellip- 
ticals split into luminosity classes. UPPER LEFT: High luminos- 
ity model ellipticals (-23.0 < Mv < -21.5); UPPER RIGHT: In- 
termediate luminosity model ellipticals (—21.5 < My < —20.5); 
LOWER LEFT: Low luminosity model ellipticals {My > -20.5); 
LOWER RIGHT: All model ellipticals. The horizontal error bar 
shows the maximum average metallicity error for model galaxies 
with sub-solar and super-solar metallicities, calculated by con- 
sidering the half-widths of the metallicity bins in the model. The 
vertical error bar shows the maximum error in the ages the model 
ellipticals. 



in age and how much by a variation in metallicity with lu- 
minosity. 

Figure |^ shows the variation in the mean ages and 
metallicities of the model cluster ellipticals with absolute 
V-band luminosity. We also show the age-metallicity param- 
eter space for these model galaxies in Figure 2] We should 
note here that the metallicity resolution in the model is low, 
with stellar mass resolved only into five metallicity bins in 
the range —1.3 < [m/H] < 0.5. We have indicated the max- 
imum average metallicity error for model galaxies with sub- 
solar and super-solar metallicities in Figure|l] by considering 



the half-widths of our metallicity bins. The resolution in age, 
by comparison, is exteremely good (0.1 Gyr), as indicated 
by the small age error bars. 

The model predicts a gradient both in the age- 
luminosity and metallicity-luminosity relations, so that 
larger ellipticals are both older and more metal-rich. 
We note first that contra ry to previous studies (e.g. 
Kauf fmann fc Charloli 1199 81) . higher mass (luminosity) el- 
lipticals are predicted to have larger mean a ges, in agree- 
ment wi th observational evid ence (see e.g. iTraeer et alJ 
l2000allbt ICaldwell et"al] |2003^ . Current understanding of 
cooling fiows in clusters is poor. Models suggest that if 
large inflows of cold gas are allowed at the centre of viri- 
alized DM haloes, it is impossible to prevent a large frac- 
tion of this material from forming stars (e.g. ICole et all 
l2000l) . which are not observed at the present epoch. To pre- 
vent this, one has to reheat or keep gas hot in massive DM 
haloes. V arious authors hav e tackl ed this problem in differ- 
ent wavs. iKauffmann et alj 1^9^, for instance, prevented 
coohng from taking place in DM haloes with circular ve- 
locities of 350 km/s and above. GALICS takes advantage 
of the observed correla tion between AGN and bulge mass 
iMaeorrian et all Il998l) and assumes that AGNs are effi- 
cient enough to completely halt cooling flows as soon as 
the bulge which harbours them reaches a critical mass of 
IO^^A/q. This coupling between AGN feedback and bulge 
mass prevents star formation early enough in large elliptical 
galaxies to allow them to grow solely through mergers of 
gas-poor progenitors. Thus, although galaxies with a larger 
mass experience their last merging events at a later time 
than their less luminous counterparts, the small gas frac- 
tion at these last-merger epochs prevents any substantial 
production of young stars from merger-driven starbursts. 
Therefore, although more massive galaxies are dynamically 
younger based on their merger record, their stellar popula- 
tions are, nevertheless, older. The average predicted age of 
a cluster elliptical is approximately 9.8 Gyrs, and the scat- 
ter in age increases towards the low mass end, in agreement 
with recent observational studies in clusters such as Virgo 
(CaldwcU ot al. 200^. The average metallicity is approxi- 
mately solar and the gradient in metallicity is modest, also 
in general agreement with r ecent spectroscopic s tudies of 
nearby cluster environments iCaldwell et alJ^2003^ . 

A c ompariso n with simple stellar population (SSP) 
models llYll2003t) shows that roughly half of the CMR 
slope is generated by the age-luminosity gradient, with the 
rest attributable to the metallicity-luminosity gradient in 
the model sample. Clearly, the age and metallicity gradi- 
ents com plement ea ch other in this model, in contrast to 
llCauffmann fc Chariot. (.1998.,) where the anti- correlation be- 
tween age and luminosity required a large compensating 
metallicity gradient (generated through high metal yields) 
to produce a CMR that was consistent with the BLE92 ob- 
servations. 

It is clearly beyond the scope of this paper to do a de- 
tailed comparison of h ow feedback is treated in GALI CS and 
the speciflc model of IKauffmann fc Chariot il998h . How- 
ever, we note that there exists at least two main differences: 



(i) Cosmological models - IKauffmann fc Chariot! (|^9^ 
adopt the SCDM cosmology while GALICS adopts the 
ACDM cosmology. 
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Figure 5. Stellar mass fraction formed at or before a given red- 
shift. The soUd line shows the cumulative mass fraction for all 
stellar mass. The other curves represent stellar mass in five dif- 
ferent metallicity ranges. 



(ii) GALICS derives feedback directly from the mass 
locked up in the spheroidal c omponent of the galaxy, while 
iKauffmann fc CharlotI il998l) use the velocity dispersion in 
dark matter haloes alone to stop the cooling. 

The first point implies that in GALICS, structures of 
a given mass will be assembled earlier on average than in 
rKauffmann fc C hariot (199^). The second point means that 
gas does not cool onto a spiral galaxy which sits in a halo 
with circular velocity greate r than or equal to 350 km/s 
in the IKauffmann fc Charlol (Il99^ model whereas it does 
in GALICS, provided the spiral does not possess a massive 
bulge. Feedback in GALICS is explicitly linked to the mass 
build up of spheroids, which in turn is correlated to the mass 
(velocity dispersion) build up of the host dark matter halo. 
However this latter correlation need not be linear, since the 
mass build up of spheroids depends on local physics (e.g. 
disk instabilitie s and mergers). We attribute the differences 
in the results of lKauffmann fc Charlo3 lll998h and GALICS 
to these two factors (cosmology and feedback modelling) but 
note that there may be other factors in the details of the 
modelling that might cause discrepancies in the relationship 
between age, metallicity and luminosity in these two models. 

Figure|5]presents the bulk cumulative SFH of the model 
cluster ellipticals. The SFH is shown both split into the five 
GALICS metallicity bins and considering all stellar mass. 
The cumulative SFH shows that 10 percent of the total 
stellar mass (solid line) was formed after z = 1, with ~ 65 
percent and ~ 40 percent already in place at z = 2 and z — 3 
respectively. The bulk SFH is quasi-monolithic because the 
low cold gas fraction at low redshifts (z < 1) means that 
merger-driven star formation does not produce substantial 
amounts of stellar material. This enables the model ellipti- 
cal CMR to maintain its slope and small scatter upto high 
redshifts (Section [IJ. 



4 EVOLUTION OF THE CMR WITH 
REDSHIFT 

We now check if it is possible to reconcile the model CMR 
with observational data at various redshifts. Figure |S] shows 
the predicted evolution of the model CMR from present day 
to a redshift of 1.27, which is roughly the redshift limit of 
current observational evidence on early-type cluster galax- 
ies (van Dokkum ct al. 2001). As before, the dotted line dis- 
plays a linear least-squares fit and we also show a progressive 
one-sigma fit to the sample, with the error bars indicating 
the local spread of points about the best-fit relation. Fig- 
ure|7|traces the evolution of the slope and the scatter. The 
shaded region denotes the area enclosed by the predicted 
slopes and their cissociated errors. 

We note that the definition of a cluster elliptical will 
change with increasing redshift. We assumed in our analy- 
sis of present-day cluster ellipticals that DM haloes with a 
mass equal to or greater than IO'^^Mq host regions of high- 
est baryonic density and therefore galaxy clusters. However, 
since DM haloes are being steadily formed through time, 
maintaining a hard mass cut-off of 10" M© for all redshifts 
would not be correct. To make this definition consistent with 
changing redshift we take into account the accretion history 
of DM haloes in the model. We first compute an average ac- 
cretion history of the present-day DM haloes with masses of 
1O^''M0 and above as a function of redshift. At each redshift 
we then define a cluster hosting DM halo as one whose mass 
is equal to or exceeds the value given by the average ac- 
cretion history. We n ote that our values are consistent with 
Ivan den BoschI ll2002l) who provides theoretical prescriptions 
for computing universal DM mass accretion histories. 

We see from Figure |7] that there is gradual evolution in 
the CMR slope, although in the range < 2 < 0.8 the evo- 
lution in the slope is zero within the errors. However, once 
we move out to z = 1.23 the change in slope is appreciable 
compared to the value at present day. Within the errors, we 
see that at high redshifts (e.g. z = 1.23) the CMR loses any 
detectable slope, partly because the expected increase in the 
scatter masks any correlation that may be present. 

In Figure |H| we put the evolution of the predicted CMR 
in the context of observational evidence. We use a variety 
of sources who have explored the CMR in various colours. 
We appl y the 30 percent co rrection for colour gradients de- 
rived bv lBower et alj lll998f) to studies that have used fixed 
apertures. We find that the slopes of the model and ob- 
served CMRs match well within the errors at all redshifts. In 
particular, we note that van Dokkum ct al. (2001) reported 
a slope at z = 1.27 that was significantly lower than the 
BLE92 value at present day. This suggests that the CMR 
slope does indeed decrease, in agreement with the expecta- 
tions of a hierarchical merger scenario. The values for the 
model scatter are also fairly consistent with the observa- 
tions, given the previous discussion in Section 2 regarding 
possible reasons for the discrepancy between the model and 
observed scatter at 2 = 0. However, we should note that the 
tightness of the predicted CMR (especially at the high lumi- 
nosity end) seems larger than what appears in observational 
studies. The model ellipticals do occupy the red part of the 
sequence (shaded region in Figure |HJ, with a scatter to bluer 
colours which increases with redshift. However, comparing 
our results at 2 ~ 0.8 to, for example Ivan Dokkum et alJ 
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Figure 6. Predicted redshift evolution of the model CMR from 
present-day to z = 1.23, which is roughly the redshift limit of 
current observational evidence on early-type cluster galaixes. Also 
shown is a linear least-squares fit (dotted line) and a progressive 
one-sigma fit, with the error bars indicating the local spread of 
points about the mean relation. 
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Figure 7. Redshift evolution of the slope and scatter in the model 
(U-V) vs. V CMR. Although the evolution in the slope is zero 
within the errors in the range < 2: < 0.8, the change in the slope 
from the value at present day becomes appreciable at z = 1.23. 
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Figure 8. CMR evolution with redshift. NOTE: We show the 
properties of CMRs in three different colours as given in the rel- 
evant studies: (U — V) data are marked as circles, (U — B) data 
are shown as triangles and (B — V) data are shown as squares. 
Filled symbols represent model values and open symbols repre- 
sent observational results. We apply the 3 percent correction for 
colour gradients derived by iBower et alj ( 1998) to studies that 
have used fixed apert ures. The obs ervational data fr om left to 
right are taken f rom: [ B ower et alj fl99^') ( marked). | S codep^p|i3 
J200ll) (markedl. Ivan Dokkum et al l 11993). lEUis et al.l ^1997^ 
Ivan Dokkum et al.l i200(l) . lvan Dokkum et alj i200ll) . We do not 
transform all results to a single colour because this requires an 
assumption of the template used to perform the transformation 
which may introduce additional uncertainities into the compari- 



('2000', Figure 8), we find that at comparable redshift, the 
observed elliptical CMR is tighter than our model predic- 
tions, although outliers do exist in the observed elliptical 
sample. 



5 PROGENITOR BIAS 

When comparing the slope and scatter of the CMR at var- 
ious redshifts, we should ideally sample the same stellar 
mass at every redshift. Only then are the slope and scat- 
ter truly meaningful tracers of the star formation history of 
the daughter mass seen today. However, an unavoidable re- 
sult of the merger paradigm is that since early-type systems 
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Figure 9. Progenitor bias: filled circles are ellipticals, triangles 
are SOs and open squares are late-type systems (spirals and irreg- 
ulars). All galaxies are progenitors of the galaxies at z = 0. The 
shaded region indicates the mean elliptical- only relation and its 
associated errors taken from Figure ISI 



form through the amalgamation of late-type units, a pro- 
gessively larger fraction of the stellar mass we see today in 
cluster ellipticals is locked up in late-type (spiral and irreg- 
ular) units at higher redshifts. Hence the early-type systems 
at high redshift form a progressively narrower subset of the 
progenitors of present-day elliptical systems. Consequently, 
by not taking into account these late-type progenitors we 
introduce a bias in the CMR, mainly in terms of the ob- 
served scatter. In this s ection we quantify the ef fect of this 
progenitor bias (see also Ivan Dokkum et al .1200 J) . Although 
tracing an astronomical object back through time is impos- 
sible observationally, it becomes a simple exercise within the 
model. 

In Figure |^ we restrict ourselves to the progenitor set 
of present-day cluster ellipticals. We show only those galax- 
ies (regardless of morphology) that eventually contribute to 
the formation of cluster ellipticals which exist at 2; = 0. We 
are therefore tracing the same stellar mass back through 
time, regardless of the type of system that hosts it. We 
find that including progenitors with SO morphology does 
not change the slope of the CMR. The overall scatter in- 
creases slightly at higher redshifts {z > 0.62), although 
the elliptical-only scatter agrees, within errors, to the E-I-SO 
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Figure 10. The number of fully formed, monolithically evolving 
elliptical progenitors of present-day cluster ellipticals at a given 
redshift. These progenitors do not undergo any further mergers, 
although quiescent star formation continues to 2 = 0. The remain- 
ing mass in the progenitor set is hosted by late- type systems. 



scatter. SOs, however, tend to contribute more outliers to the 
CMR at higher redshifts. Including the late- type progenitors 
causes the scatter to increase approximately three-fold in the 
range z > 0.8, compared to the elliptical-only scenario. This 
result agrees with CMR observations at high redshift. For 
example, van Dokkum et al. (2000) found that the elliptical 
CMR ai z — 0.83 has a scatter of ~ 0.024 while the scat- 
ter for all morphological types is ~ 0.081 - an approximate 
3.5 fold increase. Blakeslee et al. (2003) noted that for their 
observed cluster a.t z = 1.24, deriving the scatter without 
reference to morphology increases the CMR scatter three 
to four-fold. A similar increase can be estimated from the 
study by Van Dokkum et al. (2001) of a cluster ai z — 1.27 
(see their Figure 3). 

Figure [101 indicates how much of the progenitor set of 
present-day cluster ellipticals is composed of fully-formed 
ellipticals at any given redshift. It becomes clear from Fig- 
ure 1101 that if we look solely at the elliptical progenitors 
of present-day cluster ellipticals we sample a progressively 
thinner fraction of the progenitor set at higher redshift. Al- 
though restricting ourselves to this subset of progenitors 
seems to give a CMR which maintains its slope and scat- 
ter with redshift (Figure |HJ, an elliptical-only CMR can be 
used to constrain the SFH of only the part of the stellar mass 
in present-day cluster ellipticals that is contained solely in 
early-type systems at any given redshift. However, since the 
subset of elliptical progenitors is not representative of all the 
stellar material at present day, we cannot use the evolution 
of an elliptical-only CMR to constrain the SFH of the entire 
stellar mass of present-day cluster ellipticals. 

Figures |S| |5| and 1101 show that the merger paradigm 
does indeed expect to have fully formed elliptical galaxies 
(and therefore a red sequence) evolving passively at redshifts 
where CMR observations have been conducted. However, 
the elliptical-only CMR at high redshift does not correspond 
to the elliptical-only CMR at present day and comparisons 
between the two give a heavily biased picture of the star 
formation history of elliptical galaxies and has serious im- 
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plications for the ability of the CMR to discriminate between 
the monolithic collapse and hierachical merger paradigms. In 
essence, the quantity (in this case the elliptical-only CMR) 
that is being used to discriminate between the two forma- 
tion models is no longer model independent and therefore 
loses its usefulness as a discriminant. 



6 SUMMARY AND CONCLUSIONS 

We have used a semi-analytical hierarchical galaxy forma- 
tion model to investigate the existence and evolution of 
the CMR of elliptical galaxies in cluster environments. Our 
analysis shows that, by constructing a CMR purely out of 
early-type systems, the predicted relation agrees well with 
local observations (after the fixed-aperture bias has been 
corrected) and with observations at all redshifts in the range 
< z < 1.27. Secondly, we have used our analysis to quan- 
tify the issue of progenitor bias and construct the CMR that 
could be expected if we could identify all progenitor systems 
at high redshift that would eventually form part of a present- 
day cluster elliptical. We have also shown that the scatter 
in this all-progemtor CMR is consistent with the scatter 
derived, without reference to morphology, in cluster CMR 
studies at high redshift. Thirdly, we have suggested that the 
elliptical-only CMR is not a useful discriminant between the 
monolithic and merger formation scenarios since it is signif- 
icantly biased towards the monolithic picture. Although the 
merger paradigm satisfies the elliptical-only CMR in any 
case and expects to have a monolithically evolving red se- 
quence at high redshift, restricting our studies to early-type 
systems does not provide meaningful information about the 
true star formation history of all the stellar mass that is 
found today in cluster ellipticals. 

The debate regarding these two competing theories 
of elliptical galaxy formation still remains an open one. 
Although there is clear evidence of interactions, mergers 
and recent star formation in early-type systems, a possi- 
ble caveat is the inability of the merger paradigm to sat- 
isfy the high [Mg/Fe] ratios observed in luminous ellipti- 
cals fe.g. lTraeer et al.ll2000all . These super-solar abundance 
ratios indicate a lack of enrichment from Type la super- 
novae, thereby constraining the duration of star formation 
and gas infall to t imcscales shorter than a bout 1 Gyr (e.g. 
iMatteucci fc Recchi ,200L : .Ferreras fc 811^ 12003). While the 
CMR has been used as an indirect tool for constraining 
the star formation history of cluster ellipticals, more direct 
sources of evidence may be required. If the stellar mass in 
cluster ellipticals did indeed form at 2: >> 1 then we should 
not find any traces of star formation after this epoch, which 
in a ACDM universe corresponds to an age of approximately 
10 Gyrs. The merger models do of course predict star forma- 
tion right upto the present day and one could assume that at 
least a small fraction of the resultant stellar mass could be 
locked up in globular clusters, which are the faintest stellar 
aggregations that can be accessed observationally. 

Figure 1111 shows the bulk distribution of stellar mass 
in present-day cluster ellipticals predicted by the hierarchi- 
cal merger paradigm. One can treat this as a probability 
distribution of stellar mass, with the highest intensity ar- 
eas (see key) indicating ages and metallicities where most of 
the stellar mass is likely to be found. The crucial difference 




Figure 11. Predicted distribution of stellar mass contained 
in cluster ellipticals. Overplotted are o bservations o f young 
globular clusters. From left to right - [L arsen et af| j200^: 
IStrader et"an i2noaMOoudfrooii et al.N2noJl : lKissler-Patig et alJ 
j2002l) : lYi et alj 120041) . The key indicates the mass fractions cor- 
responding to the colours used in the plot. 



between this model distribution and a distribution based 
on the monolithic collapse model is the presence of young 
stars. Indeed we find that observations of young globular 
clusters ha ve be en ma de in elliptical galaxies by a variety 
of authors iGou dfrooii e t al.ll2n01 :lKissler-Patig et alJl200a : 
iLarsen et alWooa: .Strader et al.l 120031: lYi et al.ll2004^ . We 
indicate these observations in Figure lTTl The five data points 
with error bars in Figure [TTI show the age/metallicity prop- 
erties of young globular cluster populations derived in these 
studies. To conclude, we suggest that it seems increasingly 
likely that the monolithic collapse picture may simply be a 
subset of the merger paradigm and that the dominant mech- 
anism for the formation of elliptical galaxies is through the 
merging of late-type progenitors. 
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